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I. INTRODUCTION 

The Standard Model (SM) of electroweak interactions 

describes CP violation as a consequence of an irreducible 
phase in the three-family Cabibbo-Kobayashi-Maskawa 
(CKM) quark-mixing matrix [1]. In the CKM frame- 
work, tree-diagram processes dominate neutral B decays 
to CP eigenstates containing a charmonium and a if^*)*^ 
meson. These provide a direct measurement of sin2/3 [2] , 
where the angle /3 is defined in terms of the CKM matrix 
elements V^j for quarks i, j as arg[-(y^jV;*b)/(V;dV;*b)]. 

We identify (tag) the initial flavor of the reconstructed 

B candidate, -Bicc, using information from the other B 
meson, Stag) in the event. The decay rate <?+ {g-) for 
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a neutral B meson to a CP eigcnstatc / accompanied 
by a (5°) tag, before taking into account detector 
resolution effect, can be expressed as 



e-\^t\/rBO 

Atbo 



{(1tAw)±(1-2u;)x 



[Sf sin(AmdAt) - Cf cos(AmdAt)]} (1) 



where 



Cf 



2JmX 



f 



1 + IA/P' 
l-|A/f 
1 + IA/P' 



At = tree — hag is the difference between the proper de- 
cay times of Brec and -Btag, t^o is the neutral B lifetime, 
and Anid is the mass difference between the B meson 
mass eigenstates determined from B°-B° oscillations [3]. 
Here, A/ = {q/p){Af /Af) [4], where q andp are complex 
constants that relate the B-mcson flavor eigenstates to 
the mass eigenstates, and Af/Af is the ratio of the B° 
and B^ decay amplitudes to the final state /. We as- 
sume that the corresponding decay-width difference AT^ 
is zero. The average mistag probability w describes the 
effect of incorrect tags and Aw is the difference between 
the mistag probabilities for B° and B° mesons. The sine 
term in Eq. 1 results from the interference between di- 
rect decay and decay after B^-B^ oscillation. A non-zero 
cosine term arises from the interference between decay 
amplitudes with different weak and strong phases (direct 
CP violation ^ 1) or from CP violation in 

B^ mixing {\q/p\ ^ 1). In the SM, CP violation in mixing 
and direct CP violation are both negligible in 6 — > ccs de- 
cays [4]. Under these assumptions, Aj = rjfe~'^^^ , where 
r]f = +\ (—1) is the CP eigenvalue for a CP-even (odd) 
final state, implying C/ = 0. Thus, the time-dependent 
CP-violating asymmetry is 



Acp{At) 



_ g+{At) - g.{At) 
g+iAt) + g-iAt) 
= {l~2w)SfSm{AmdAt), 



(2) 



and Sf = — 77/sin2/3. If the assumption that C/ = is 



relaxed, then S r = — r 



1 



C? sin2/3. 



'f = -Vfy---^f' 
In a previous publication [5], we reported time- 
dependent CP asymmetries in terms of the parameters 
sin2/3 and |A/|. In this paper we report results in terms 
of Sf and C/ to be consistent with other time-dependent 
CP asymmetry measiirements. We reconstruct B" de- 
cays to the final states J/^K°, J/^pK^, ipi2S)K°, XciK°, 
r]cK% and J/ijK*{Smf with K*{Smf K°7rO [6]. The 
J/ip Kl final state is CP-even and the J/ip K*'^ final state 
is an admixture of CP-even and CP-odd amplitudes. The 
remaining final states are CP-odd. The CP-even and odd 
amplitudes in i?" J/ipK*^ decays can be separated in 
an angular analysis [7] . In this analysis, we average over 



the angular information resulting in a dilution of the mea- 
sured CP asymmetry by a factor 1 — 2_Rj^, where is 
the fraction of the L=l contribution. In Rcf. [7] we have 
measured i?_L = 0.233 ± 0.010 (stat) ± 0.005 (syst), which 
gives an effective r?/ = 0.504 ± 0.033 after acceptance 
corrections for / = J/tpK*'^. In addition to measuring a 
combined Sf and Cf for the CP modes described above, 
we measure Sf and Cf for each final state / individu- 
ally. We split the J/ipK^ mode into samples with cither 
ifg ^ TT+TT" or 7r°7r°. We also combined J/tpK'^ channel 
with either a Kg or K^. Compared to our previous 
publication [5], the current analysis contains 82 x 10^ ad- 
ditional BB decays and improved track reconstruction 
algorithms have been applied to the entire data set. 



II. THE DATA SET AND BABAR DETECTOR 

The results presented in this paper are based on data 
collected with the BaBar detector at the PEP-II asym- 
metric energy storage rings [8] operating at the 
SLAC National Accelerator Laboratory. At PEP-II, 9.0 
GeV electrons and 3.1 GeV positrons collide at a center- 
of-mass energy of 10.58 GeV, which corresponds to the 
mass of the Y(AS) resonance. The asymmetric energies 
result in a boost from the center-of-mass (CM) frame to 
the laboratory of (i-y « 0.56. The data set analyzed has 
an integrated luminosity of 425.7 fb~^ corresponding to 
(465± 5) X 10® BB pairs recorded at the T(4S') resonance. 

The BaBar detector is described in detail elsewhere [9] . 
Surrounding the interaction point is a five-layer, double- 
sided silicon vertex tracker (SVT), which measures the 
impact parameters of charged particle tracks in both 
the plane transverse to, and along, the beam direction. 
A 40-layer drift chamber surrounds the SVT and pro- 
vides measurements of the momenta for charged parti- 
cles. Charged hadron identification is achieved through 
measurements of particle energy-loss in the tracking sys- 
tem and the Cherenkov angle obtained from a detector 
of internally reflected Cherenkov light. A CsI(Tl) elec- 
tromagnetic calorimeter (EMC) provides photon detec- 
tion, electron identification, and reconstruction. The 
aforementioned components arc enclosed by a solenoid 
magnet, which provides a 1.5 T magnetic field. Finally, 
the flux return of the magnet (IFR) is instrumented in 
order to allow discrimination of muons from pions. For 
the most recent 211. 7fb^^ of data, a portion of the re- 
sistive plate chambers in the IFR has been replaced by 
limited streamer tubes [10]. 

We use a right-handed coordinate system with the z 
axis along the electron beam direction and the y axis up- 
ward. Unless otherwise stated, kinematic quantities are 
calculated in the laboratory rest frame. We use Monte 
Carlo (MC) simulated events generated with the BaBAR 
simulation based on GEANT4 [11] for detector responses 
and EvtGen [12] for event kinematics to determine signal 
and background characteristics, optimize selection crite- 
ria, and evaluate efficiencies. 
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III. RECONSTRUCTION OF B CANDIDATES 

We select two samples of events in order to measure 

the timc-dcpcndcnt CP asymmetry parameters Sf and 
Cf. a sample of signal events used in the extraction of 
the CP parameters (Bcp) and a sample of fully recon- 
structed B meson decays to flavor cigcnstates (Bflav)- 
The Bcp sample consists of B*^ decays to J/tpK^, J/ipK^, 
V'(2S')if°, ricK^, XciK^ and J/^pK*{892f, where K*° 
decays to K'^tt'^. The i?flav sample consists of i?° decays 
to _D(*'^(7r+, p+, a^) final states. We use the i?flav sam- 
ple to determine the dilution (mistag probability) and 
the resolution function, discussed in Section V. We as- 
sume that the interference between the CP side and the 
tag side reconstruction is negligible and therefore that 
the dilution and resolution parameters are the same for 
the -Bflav and Bcp samples. We also select a sample of 
fully reconstructed charged B meson decays to J/tpK~^, 
i;{2S)K+, XciK+, rjcK+, and J/i)K*{892)+ , where K*+ 
decays to K~^n^ or K^tt'^, to use as a control sample. 

The event selection is unchanged from that described 
in Ref [5]. J/tp and tp{2S) mesons are reconstructed via 
their decays to e~^e~ or final states. At least one of 

the leptons is required to pass a likelihood particle identi- 
fication algorithm based on the information provided by 
the EMC, the IFR and from ionization energy loss mea- 
sured in the tracking system. We require the invariant 
mass of the muon pair m(/Lt+/Lt~) to be in the mass range 
3.06 3.14 GeV/c2 for J/tP or 3.636-3.736 GeV/c^ for iP{2S) 
candidates. For J/tp e+e" and tp{2S) e+e" decays, 
where the electron may have radiated bremsstrahlung 
photons, part of the missing energy is recovered by iden- 
tifying neutral clusters with more than 30 MeV lying 
within 35 mrad in polar angle and 50 mrad in azimuth 
of the electron direction projected onto EMC. The in- 
variant mass of e+e~ pairs is required to be within 2.95- 
3.14 GeV/c2 for J/^P candidates, or 3.436-3.736 GeV/c^ for 
ip{2S) candidates. 

We also construct ip{2S) mesons in the J/'^tt+tt" final 
state, where the J/ip candidate is combined with a pair 
of oppositely-charged tracks assumed as pions with no 
particle identification applied, and the pion pair-invariant 
mass between 400MeV/c2 and 600 MeV/ c^. Candidates 
with 3.671 GeV/c2 < m{J/tpTT+TT-) < 3.701 GeV/c^ are 
retained. 

The Xci candidates are reconstructed in the J/tp J final 
state. The photon candidates are required to have an 
energy greater than 100 MeV but less than 2 GeV, and, 
when combined with other photons, not to form a 7r° 
candidate with invariant mass 120 MeV/c^ < m{j'y) < 
150 MeV/c^. The; invariant mass of the Xci candidate is 
required to be between 3.477 GeV/c^ and 3.577 GeV/c^. 
Mass constraints are applied in the fits to improve the 
determinations of the energies and momenta of the J/tp , 
tp{2S), and Xci candidates. 

We reconstruct the B^ VcK^ mode using the 
rjc — > K'^K^tt~ decay. We exploit the fact that the 
r]c decays predominantly through a Kn resonance at 



around 1.43 GeV/c^ and a K';!K resonance close to the 
threshold. We require that 'm{K'^7r^) or •m{K^iT~) is 
within mass range of 1.26 GcV/c^ and 1.63 GeV/c^, or 
1.0 GeV/c2 < m{K+K^s) < 1.4 GeV/c^. 

The decay channels if+7r^, K^tt^tt''^ if "'"7r"'"7r~7r~, 
and i^gTr+TT" are used to reconstruct , while D~ can- 
didates are selected in the K~^'K~n~ and Kgn~ modes. 
We require that the and D~ candidate invariant 
mass is within ±3(7 of their respective nominal mass, 
where a is the uncertainty calculated for each candidate. 
A mass-constrained fit is then applied to the D'^ and 
D~ candidates satisfying these requirements. We form 
D*~ candidates in the decay D*^ D°7r~ by com- 
bining a with a pion that has momentum greater 
than 70 MeV/c. The D*~ candidates are required to have 
m{D'^TT^) within ±1.1 MeV/c^ of the nominal D*~ mass 
for the D° K+-k--k^ mode and ±0.8 MeV/c^ for all 
other modes. 

For the J/tpK^ decay, we use both K° tt+tt" and 
decays; for other B decay modes we only 
use K° TT+TT". Candidates in the Kg tt+tt" mode 
are selected by requiring an invariant 7r"'"7r~ mass, com- 
puted at the vertex of the two oppositely-charged tracks, 
between 472.67 MeV/c^ and 522.67 MeV/c^. We further 
apply a mass constraint fit to the K'^ candidates before 
combining them with charmonium candidates to form 
B^ candidates. Neutral pion candidates, in the mass 
range 100-155 MeV/c^, are formed from two 7 candidates 
from the EMC. Pairs of 7r° are combined to construct 
—> tt'^tt'^ candidates. The minimum energy is required 
to be 30 MeV for 7, 200 MeV for n°, and 800 MeV for K° 
candidates. To select K'^ candidates, the 7r°7r° invariant 
mass is restricted to the region between 470 MeV/ and 
550MeV/c2. 

Candidates for are identified in the EMC and IFR 
detectors as reconstructed clusters that cannot be associ- 
ated with any charged track in the event. As the energy 
of cannot be measured well, the laboratory momen- 
tum of the K'l is determined by its flight direction and the 
constraint that the invariant mass of the J/tpK^ system 
has the known B^ mass. For events with multiple J/tpK^ 
candidates, a hierarchy is imposed where the highest en- 
ergy EMC cluster for multiple EMC combinations, or the 
IFR cluster with the largest number of layers for multi- 
ple IFR combinations, is selected. In case both EMC 
and IFR combinations are found, the EMC combination 
is chosen because of its better angular resolution. 

We reconstruct K*'^ candidates in the i^sTr*^ mode, 
while K*'^ candidates are reconstructed in the if +7r" and 
if^7r+ modes. The invariant mass of the two daughters 
is required to be within ±100 MeV/c^ of the nominal K* 
mass. 

The /9+ candidates are reconstructed in the tt+tt^ fi- 
nal state, where the tt+tt*^ mass is required to lie within 
±150 MeV/c^ of the nominal p"*" mass. Candidates in the 
decay mode — > 7r+7r~7r+ are reconstructed by combin- 
ing three charged tracks with pion mass assumption, and 
restricting the three-pion invariant mass to lie between 
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1.0 and 1.6GcV/c2. 

Events that pass the selection requirements are re- 
fined using kinematic variables. For the J/ipKl mode, 
the difference /S.E between the candidate CM energy 
and the beam energy in the CM frame, i^beann is re- 
quired to satisfy \^E\ < 80 MeV. For all other categories 
of events, we require \AE\ < 20 MeV and the beam- 
energy substituted mass rriEs = v^(-^bcam)^ - (pW^ to 
be greater than 5.2GeV/c^, where p*g is the B momen- 
tum in the CM frame. When multiple B candidates (with 
Toes > 5.2 GeV/c^) are found in the same event, the can- 
didate with the smallest value of \AE\ is selected. 

We calculate the proper time difference At between the 
two B decays from the measured separation Az between 
the decay vertices of B^cc and Stag along the collision 
(z) axis [13]. The z position of the i?roc vertex is de- 
termined from the charged daughter tracks. The Stag 
decay vertex is determined by fitting tracks not belong- 
ing to the Sroc candidate to a common vertex, including 
constraints from the beam spot location and the B^^c 
momentum [13]. Events are accepted if the calculated 
At uncertainty is less than 2.5 ps and \At\ is less than 
20 ps. The fraction of signal MC events satisfying such a 
requirement is 95 %. 

IV. B MESON FLAVOR TAGGING 

A key ingredient in the measurement of time- 
dependent CP asymmetries is the determination of 
whether the Brec was a or a B^ at the time of At = 0. 
This 'flavor tagging' is achieved with the analysis of the 
decay products of the recoiling B meson Stag- The over- 
whelming majority of B mesons decay to a final state that 
is flavor-specific, i.e., only accessible from either a B^ or 
a B^. The purpose of the flavor-tagging algorithm is to 
determine the flavor of Stag with the highest efficiency 
etag and lowest probability w of assigning the wrong fia- 
vor. It is not necessary to fully reconstruct Stag in order 
to determine its fiavor. 

The figure of merit for the performance of the tagging 
algorithm is the effective tagging efficiency 

Q = etag(l - 2wf, (3) 

which is related to the statistical uncertainty as and ac 
in the coefficients Sf and C/ through 

as,c oc (4) 

The tagging algorithm we employ [5, 13] analyzes 
tracks on the tag side to assign a flavor and associated 
probability to -Stag- The flavor of Stag is determined from 
a combination of nine different tag signatures, such as iso- 
lated primary leptons, kaons and pions from B decays to 
final states containing D* mesons, and high momentum 
charged particles from B decays. The properties of those 
signatures are used as inputs to a single neural network 



that is trained to assign the correct flavor to Stag- The 
output of this neural network then is divided into seven 
mutually-exclusive categories. These are (in order of de- 
creasing signal purity) Lepton, Kaon I, Kaon II, Kaon- 
Pion, Pion, Other and Notag. The events with the neural 
network output \NN\ > 0.8 are defined as Lepton cate- 
gory, if they are also accompanied by an isolated primary 
lepton; otherwise they are categorized as Kaon / tag. For 
the other five tag categories (Kaon II, KaonPion, Pion, 
Other and Notag) the outputs of the neutral network are 
required to satisfy: 0.6 < |iViV| < 0.8, 0.4 < \NN\ < 0.6, 
0.2 < \NN\ < 0.4, 0.1 < \NN\ < 0.2, and \NN\ < 0.1, 
respectively. 

The performance of this algorithm is evaluated using 
the -Bflav sample. The final state of the -Bflav sample can 
be classified as mixed or unmixed depending on whether 
the reconstructed flavor-eigenstate Sfiav has the same or 
opposite fiavor as the tagging B. After taking the mistag 
probability into account, the decay rate g±^B" {g± go) 
for a neutral B meson to decay to a flavor eigenstate 
accompanied by a B° (5°) tag can be expressed as 

g±,Bo{At) a [{1 - Awi) ± {1 - 2wi) cos{AmdAt)] 

(5) 

g^^Boi^t) oc [(l + Aw,)±(l-2w,)cos(AmdAi)], 

where the ± sign in the index refers to mixed (— ) and 
unmixed (-I-) events; the index i denotes the ith tagging 
catcgoiy. The performance of the tagging algorithm is 
summarized in Table I. The events in the Notag cate- 
gory contain no flavor information, so carry no weight 
in the time-dependent analysis. They are excluded from 
further analysis. The total effective tagging efficiency is 
measured to be (31.2 ± 0.3)%. 

TABLE I: Efficiencies ei, average mistag fractions Wi, mistag 
fraction differences between and -S" tagged events Awi, 
and effective tagging efficiency Qi extracted for each tagging 
category i from the Bflav sample. 



Category ei (%) Wi (%) Awi (%) Qi (%) 

Lepton 8.96 ±0.07 2.8 ± 0.3 0.3 ± 0.5 7.98 ±0.11 

Kaon I 10.82 ± 0.07 5.3 ± 0.3 -0.1 ± 0.6 8.65 ± 0.14 

Kaon 11 17.19 ±0.09 14.5 ±0.3 0.4 ±0.6 8.68 ± 0.17 

KaonPion 13.67 ± 0.08 23.3 ± 0.4 -0.7 ± 0.7 3.91 ± 0.12 

Pion 14.18 ±0.08 32.5 ±0.4 5.1 ±0.7 1.73 ±0.09 

Other 9.54 ±0.07 41.5 ±0.5 3.8 ± 0.8 0.27 ±0.04 

All 74.37 ± 0.10 31.2 ± 0.3 



V. LIKELIHOOD FIT METHOD 

We determine the composition of our final sample by 

performing simultaneous fits to the m^g distributions for 
the full Bcp and Sfiav samples, except for the J/tpK^ 
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sample for which wc extract the K^^ momentum by using 
the mass constraint and fit the /S.E distribution. Wc 
then perform a simultaneous maximum likelihood fit to 
the At distribution of the tagged Bcp and .Bflav samples 
to measure Sf and C/. 

We define a signal region of 5.27 < toes < 5.29 GeV/c^ 
{\IS.E\ < lOMeV for JjiiKl), which contains 15481 can- 
didate events of Bcp sample that satisfy the tagging and 
vcrtcxing requirements (sec Table II). The signal toes 
distribution for the full Bqp and Sflav samples, except for 
the J/xjjK^ sample is described by a Gaussian function. 
The background toes distribution is modeled by an AR- 
GUS threshold function [14], where a shape parameter is 
allowed to vary in the fit. For the decay modes of J/iIj if", 
ijj2sKg, XciKg, J/iIjK*° and Bflav, we use simulated 
events to estimate the fractions of background events that 
peak in the toes signal region (toes > 5.27 GeV/c^) due 
to cross- feed from other decay modes. We describe this 
component with a Gaussian function having the same 
mean and width as the signal and refer to it as peak- 
ing background because if neglected, it would lead to 
an overestimate of the signal yield. The peaking back- 
ground is less than 1 % in the decay of J/tpK'^, and 
at the level of a few percent in most other decay modes. 
The only exception is the decay of i?" J/tpK*'^, where 
the peaking background level is about 13 %. MC simula- 
tions show that it consists of 44 % of i?"*" decays, 32 % of 
B° XcK° decays and 24 % of other B" decays. For the 
ricK'^ mode, the cross-feed fraction is determined from a 
fit to the rriKKTr and toes distributions in data. For the 
J/ijjKl decay mode, the signal AE distribution is deter- 
mined from MC simulated events. The sample composi- 
tion, effective r]f, and AE distribution of the individual 
background sources are determined either from simula- 
tion (for B — !■ J/ipX) or from the mg+g- sidebands in 
data (for non-J/i/i background). Figure 1 shows the dis- 
tributions of Toes obtained for the Bcp and Bflav events, 
and AE obtained for the J/i/jK^ events. 

The At distributions of the Bcp sample are modeled by 
Eq. 1 and those of the Bflav sample by Eq. 5. The At dis- 
tributions for the signal are convolved with a resolution 
function common to both the Bflav and Bcp samples, 
modeled by the sum of three Gaussian functions [13], 
called the core, tail and outlier components, which can 
be represented as a function of the reconstruction uncer- 
tainty St = At — Attruc as follows: 

'^coroO'At, S'corcO'At) 
+ /tail/lG(<5<; ^tail^Af, StailCTAf) 
+ /out/lG(^i;^out,5'out), (6) 



where 



and 



(8) 



/core ~^ /tail ~t~ /out — f- 



The widths (cr) of the core and tail components include 
two independent scale factors, S'corc and 5taii, to accom- 
modate an overall underestimate or overestimate of the 
At measurement error a At for all events. The parame- 
ter Score is free in the fit and its value is close to unity. 
The value of S'taii is derived from MC studies and fixed 
to be 3. Studies show that the measurement of C/ and 
Sf is not sensitive to the choice of the 5taii value. We 
later vary the S'taii value within a large region and as- 
sign the shift of the measured C/ and Sf values as the 
corresponding systematic uncertainties. We account for 
residual charm decay products included in the Btag can- 
didate vertex by allowing the core and tail Gaussian func- 
tions to have non-zero mean values (bias. Score and 
<5taii 7^ 0). The bias (5core) and width (Score) of the core 
component are allowed to differ for the lepton-tagged and 
nonlepton-tagged events. We use common parameters for 
the tail component. In order to account for the strong 
correlations with other resolution parameters, the outlier 
component bias (5out) and width (Sout) are fixed to ps 
and 8 ps, respectively. 

The At spectrum of the combinatorial background is 
described by an empirical distribution, consisting of com- 
ponents with zero and non-zero lifetimes (rbg) that are 
convolved with a resolution function [13] distinct from 
that used for the signal. Here, we use a double-Gaussian 
distribution, which has components similar to the core 
and outlier distributions described above. In this case, 
the resolution function is common to all tagging cate- 
gories. The peaking background is assigned the same At 
distribution as the signal but with Sf = Cf = 0, and uses 
the same At resolution function as the signal. The non- 
zero lifetime component of the combinatorial background 
contains both mixed and un- mixed events. Therefore we 
allow the value of Ato^ for this component {Amd,hg) to 
vary in the fit. 

In addition to Sf and Cf, there are 69 free parameters 
in the fit. For the signal, these are 

• 7 parameters for the At resolution: i5coro and Score 
for the lepton-tagged and nonlepton-tagged events, 

/core) /tail) and 

^tail- 

• 12 parameters for the average mistag fractions 

and the differences Awi between B° and B^ mistag 
fractions for each tagging category, 

• 1 parameter for the small difference between B° 
and B^ reconstruction efficiency [13], and 

• 6 parameters for the small difference between B° 
and B° tagging efficiencies in each tagging cate- 
gory [13]. 

The background parameters that are allowed to vary are 

• 24 mistag fraction parameters: Wi and Awi of each 

tagging category for background components with 
zero and non-zero lifetime, respectively. 
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FIG. 1: Distributions for Bcp and Bflav candidates satisfying the tagging and vertexing requirements: a) mES for the final 
states J/ijK°, ij{2S)K°, XaK'^s, and ricK°, b) AE for the final state J/^jKI, c) tues for J/iP K'" (K*'^ K°-iv°), and d) rriES 
for the Bflav sample. In each plot, the shaded region is the estimated background contribution. 



• 3 parameters for the At resolution: ^corei 'S'core 

/core* 



and 



4 parameters for the Bflav time dependence: 2 pa- 



rameters for the fraction (/j 



prompt. 



of zero lifetime 



component for the lepton-tagged and nonlepton- 
tagged events, rbg, and Amd^bg- 

• 8 parameters for possible CP violation in the back- 
ground, including the apparent CP asymmetry of 
non-peaking events in each tagging category, 

• 1 parameter for possible direct CP violation in the 
XciK^l background coming from J/ifjK*^ , and 

• 3 parameters for possible direct CP violation in the 



J/tjjKl mode, coming from J/tpK^, J/tpK*^, 
the remaining J/il) backgrounds. 



and 



The effective value of |A/| for the non-J/ijj background is 
fixed from a fit to the J/V' -candidate sidebands in Jf-ipK^. 
We fix tb« = 1.530 ps and Anid = 0.507 ps"! [3]. The 
determination of the mistag fractions and At resolution 
function parameters for the signal is dominated by the 
Bfiav sample, which is about 10 times larger than the CP 
sample. 
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VI. LIKELIHOOD FIT VALIDATION 

We perform three tests to validate the fitting proce- 
dure. The first of these tests consists of generating en- 
sembles of simulated experiments from the probability 
density function (PDF) and fitting each simulated exper- 
iment. We determine that the fitted values of Sf and C/ 
parameters are unbiased, and that the fit returns reason- 
able estimates of the statistical uncertainties, by verifying 
the distribution of the pull P on a parameter O, given by 
V = (Ofit — C'gcn)/o'(Cfit), is consistent with a Gaussian 
centered about zero with a width of one. The quantity 
Ofit is the fitted value, with a fitted error of a{Ofit), and 
Ogen is the generated value. 

The second test involves fitting simulated signal events 
that include the full BaBAR detector simulation. For each 
decay mode, we divide the signal MC sample to many 
data-sized samples, fit them one by one, and then exam- 
ine the distribution of the fitted results. We make sure 
that the V distributions for these signal-only simulated 
experiments are consistent with a Gaussian distribution 
centered at zero with a width of one. 

The third test is to perform null tests on control sam- 
ples of neutral and charged B events where Sf and Cf 
should be very small or zero. The parameters Sf and 
Cf are consistent with zero for the charged B sample 
of J/^K^, ^{2S)K^, XciK^, and J/ipK*^ final states. 
For the neutral i?flav sample, we find that the Sf and 
Cf parameters slightly deviate from zero at approxi- 
mately twice the statistical uncertainty (see Table II). 
The deviation of 5*/ from zero is consistent with the 
directly measured CP asymmetry S ^ — 2rsin(2/3 -|- 
7)cos(^) < 0.04 [15] in 5° i:»(*)±/i=F [16] due to in- 
terference from doubly-CKM-suppressed decays, where 
7 = arg[-(KidKrb)/(^cdKb)]' ^ is the strong phase differ- 
ence between CKM-favored and doubly-CKM-suppressed 
amplitudes, and r ~ 0.02 is the ratio of the two am- 
plitudes. Considering this expected CP asymmetry in 
the -Bfiav sample and systematic uncertainties (at ~ 1% 
level), we conclude that our analysis is free of pathologi- 
cal behaviors. 



VII. RESULTS 

The fit to the Bcp and -Bflav samples yields —T]fSf = 
0.687 ± 0.028 and Cf = 0.024 ± 0.020, where the errors 
are statistical only. The correlation between these two 
parameters is -1-0.1 %. We also performed the fit using 
sin2/3 and |A/| as fitted parameters, and found sin2/3 — 
0.687 ± 0.028 and |A/| = 0.977 ± 0.020. The correlation 
between the fitted sin2/3 and |A/| parameters is —0.14%. 
Figure 2 shows the At distributions and asymmetries in 
yields between events with B^ and B° tags for the ??/ = 
— 1 and rjf = +1 samples as a function of At, overlaid 
with the projection of the likelihood fit result. Figure 3 
shows the time-dependent asymmetry between unmixed 
and mixed events for hadronic B candidates with toes > 




FIG. 2: a) Number of ry/ = -1 candidates {J/ipKl, i,{2S)K°s, 
XciKs, and rjcK'^) in the signal region with a B° tag (Nga) 
and with a B° tag (Ngo), and b) the raw asymmetry, 
(Ngo — Ngo) /{Ngo + Ngo), as functions of At; c) and d) 
are the corresponding distributions for the rjf — +1 mode 
Jf-tpKl. The solid (dashed) curves in (a) and (c) represent 
the fit projections in At for B'^ (B") tags. The shaded regions 
represent the estimated background contributions to (a) and 
(c). The curves in (b) and (d) are the fit projections of the 
raw asymmetry between tagged and B^ tagged events. 



5.27 CeV/c^. We also perform a fit in which we allow 
different 5*/ and Cf values for each charnionium decay 
mode, a fit to the J/'ipK'^ (tt+tt" -t- tt^tt^) mode, and a fit 
to the J/iIjK^ {K° + Kl) sample. The results for some 
of these studies are shown in Figure 4. We split the 
data sample by run period and by tagging category. We 
perform the CP measurements on control samples with 
no expected CP asymmetry. The results of these fits are 
summarized in Table II. 

The dominant systematic uncertainties on Sf are sum- 
marized in Tables III and IV. The dilution due to flavor 
tagging can be different between Bcp and -Bflav events. 
We study this effect by comparing the results in large 
samples of simulated Bcp and -Bflav events. The uncer- 
tainties due to At resolution functions for both signal 
and background components are estimated by varying 
the fixed parameters and by using alternative models. We 
also vary the peaking background fractions based on esti- 
mates derived from simulation, and vary the CP content 
of the background over a wide range to estimate the effect 
due to our limited knowledge of background properties. 
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TABLE II: Number of events A'tag and signal purity P in the signal region after tagging and vertexing requirements, and results 
of fitting for CP asymmetries in the Bcp sample and various subsamples. Fit results for the Bflav and control samples are 
also shown here. Errors are statistical only. 



Sample 


Ntag 


P(%) 


—rifSf 

U J 


J 


FuU CP sample 


15481 


76 


0.687 ± 0.028 


0.024 ± 0.020 


J/lllKli-K+TV') 


5426 


96 


0.662 ±0.039 


0.017 ± 0.028 




1324 


87 


0.625 ±0.091 


0.091 ± 0.063 




861 


87 


0.897 ±0.100 


0.089 ± 0.076 




385 


88 


0.614 ±0.160 


0.129 ±0.109 




381 


79 


0.925 ±0.160 


0.080 ±0.124 




5813 


56 


0.694 ±0.061 


-0.033 ± 0.050 




1291 


67 


0.601 ± 0.239 


0.025 ± 0.083 




6750 


95 


0.657 ± 0.036 


0.026 ± 0.025 




12563 


77 


0.666 ± 0.031 


0.016 ±0.023 


77/ = —1 


00 ( ( 


yo 






1999-2002 data 


3079 


78 


0.732 ±0.061 


0.020 ± 0.045 


2003-2004 data 


4916 


77 


0.720 ±0.050 


0.045 ± 0.036 


2005-200fi data 

^\J\J\J ^\J\J\J l-lChKiCb 


4721 


76 


632 ± 052 


027 ± 037 


2007 data 


2765 


75 


0.663 ±0.071 


-0.023 ± 0.049 


Lepton 


1740 


83 


0.732 ± 0.052 


0.074 ± 0.038 


Kaon I 


2187 


78 


0.615 ±0.053 


-0.046 ± 0.039 


Kaon II 


3630 


76 


0.688 ±0.056 


0.068 ± 0.039 


KaonPion 


2882 


74 


0.741 ± 0.086 


0.013 ±0.061 


Pion 


3053 


76 


0.711 ±0.132 


0.016 ± 0.090 


Other 


1989 


74 


0.766 ± 0.347 


-0.176 ±0.236 


-Bflav sample 


166276 


83 


0.021 ±0.009 


0.012 ±0.006 


B"*" sample 


36082 


94 


0.021 ±0.016 


0.013 ±0.011 




At(ps) 



FIG. 3: Time-dependent asymmetry between unmixed and 
mixed events, A^^^ = (A^unmix " ^mix)/(^unmix + ^mix). 
as a function of At for hadronic B candidates (Bflav) with 
TUES > 5.27 GeV/c^ The curve is the corresponding fit pro- 
jection. 



The uncertainties in the J/ipKl sample are studied by 
varying the compositions of the signal and background, 
by modifying the AE PDF based on studies performed 
with the J/ipK'^ control sample, and by varying the 
branching fractions of the background modes and their 
CP asymmetries. Other sources of uncertainty such as 



the values of the physics parameters Am^, tb, AFd/r^, 
the beam spot and detector alignment, and other fixed 
parameters, are studied by varying them according to 
their world averages, the calibration, and the statistical 
uncertainty, respectively. Despite the large amount of 
simulated signal events that included the full BaBAR de- 
tector simulation, we can only validate the possible fit 
bias to be no more than certain precision. As a result, 
we assign a systematic uncertainty corresponding to any 
deviations and the statistical uncertainties of the mean 
values of the fitted Sf and C/ from the generated values 
as the possible fit bias (MC statistics). 

The only sizable systematic uncertainties on C/ are 
due to the CP content of the peaking backgrounds and 
due to the possible interference between the suppressed 
b — > ucd amplitude with the favored b cud amplitude 
for some tag-side B decays [15]. The total systematic 
error on Sf (C/) is calculated by adding the individual 
systematic uncertainties in quadrature and is found to 
be 0.012 (0.016). The main sources of systematic uncer- 
tainty are listed in Tables III and IV. 

For the rjcK^ mode, we found —rjfSf — 0.925 ± 
0.160(stat)±0.057(syst), which has a significance of 5.4 ct 
standard deviations including systematic uncertainties. 
Our result is the first observation of CP violation in this 
mode. 
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At (ps) 



At (ps) 



FIG. 4: a) Number of candidates in the signal region with a 5° tag {Ngo ) and with a tag (Ngo), and b) the corresponding 
raw asymmetry, {Ngo — Ngo)/{Ngo + Ngo), as functions of At for each B'^ decay modes. The solid (dashed) curves represent 
the fit projections in At for _B" (B") tags. The right-hatched (left-hatched) shaded regions represent the estimated background 
contributions in (B") tagged events. 



VIII. CONCLUSIONS 



find 



We report improved measurements of the time- 
dependent CP asymmetry parameters. The results in 
this paper supercede those of our previous publication [5] . 
We report our measurements in terms of C/ and Sf. We 



Cf = 0.024 ± 0.020 (stat)± 0.016 (syst), 
-r]fSf = 0.687 ± 0.028 (stat)± 0.012 (syst), 

providing an independent constraint on the position of 
the apex of the Unitarity Triangle [17]. Our mea- 
surements agree with previous published results [5, 18] 
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TABLE III: Main systematic uncertainties on Sf and C/ for the full CP sample, and for the J/ip K'^ , J/tpKs, and J/ip K° 
samples. For each source of systematic uncertainty, the first line gives the error on Sf and the second line the error on C/. 
The total systematic error (last row) also includes smaller effects not explicitly mentioned in the table. 



Source/sample 




Full 


J/%I}K° 


J/^Kl 




Beam spot 


Sf 


0.001 


0.002 


0.003 


0.000 




Cf 


0.001 


0.001 


0.002 


0.000 


Mistag differences 


Sf 


0.006 


0.006 


0.006 


0.006 




Cf 


0.002 


0.002 


0.002 


0.002 


At rcsohition 


Sf 


0.007 


0.007 


0.007 


0.007 




Cj 


().()():i 


().()();', 




0.007 


J/ip K'^ background 


Sf 


0.006 


0.006 


0.000 


0.027 




Cf 


0.001 


0.001 


0.000 


0.004 


Backgrouiid fraction 


Sf 


0.005 


0.004 


0.004 


0.004 


and CP content 


Cf 


0.003 


0.002 


0.001 


0.011 


rriES parameterization 


Sf 


0.002 


0.002 


0.003 


0.001 




L'/ 


U.UUU 


U.UUi 


U.UUi 


n nnn 
U.UUU 


Amd,rs, Ard/Fd 


Sf 


0.003 


0.003 


0.004 


0.004 




Cf 


0.001 


0.001 


0.001 


0.001 


Tag-side interference 


Sf 


0.001 


0.001 


0.001 


0.001 




Cf 


0.014 


0.014 


0.014 


0.014 


Fit bias 


Sf 


0.002 


0.004 


0.004 


0.006 


(MC statistics) 


Cf 


0.003 


0.004 


0.004 


0.006 


Total 


Sf 


0.012 


0.013 


0.012 


0.031 




Cf 


0.016 


0.018 


0.016 


0.027 



TABLE IV: Main systematic uncertainties on Sf and C/ for the J/i/) A'2(7r+7r~), J/i/) A'g(7r°7r°), '>1}{2S)K%, Xci^^^s, n^K%, and 
J/'il)K*'^{K*^ — > ifg""") decay modes. For each source of systematic uncertainty, the first line gives the error on 5*/ and the 
second line the error on Cf. The total systematic error (last row) also includes smaller effects not explicitly mentioned in the 
table. 



Source/sample 




J/tpK°{TV+-K- 




^i2S)K° 






J/i!K*° 


Beamspot 


Sf 


0.003 


0.002 


0.008 


0.028 


0.001 


0.006 




Cf 


0.002 


0.003 


0.009 


0.012 


0.000 


0.000 


Mistag differences 


Sf 


0.006 


0.006 


0.006 


0.006 


0.006 


0.006 




Cf 


0.002 


0.002 


0.002 


0.002 


0.002 


0.002 


At resolution 


Sf 


0.007 


0.007 


0.007 


0.010 


0.016 


0.026 




Cf 


0.003 


0.004 


0.007 


0.004 


0.004 


0.006 


J/tpKl background 


Sf 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 




Cf 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


Background fraction 


Sf 


0.003 


0.007 


0.016 


0.017 


0.051 


0.056 


and CP content 


Cf 


0.001 


0.003 


0.006 


0.010 


0.019 


0.026 


mES 


Sf 


0.002 


0.009 


0.024 


0.006 


0.002 


0.037 


parameterization 


Cf 


0.001 


0.006 


0.001 


0.002 


0.001 


0.008 


Amd,rs,Ard/rd 


Sf 


0.003 


0.007 


0.016 


0.003 


0.016 


0.014 




Cf 


0.001 


0.001 


0.001 


0.001 


0.002 


0.001 


Tag-side interference 5/ 


0.001 


0.001 


0.001 


0.001 


0.001 


0.001 




Cf 


0.014 


0.014 


0.014 


0.014 


0.014 


0.014 


Fit bias 


Sf 


0.005 


0.004 


0.008 


0.007 


0.007 


0.027 


(MC statistics) 


Cf 


0.004 


0.003 


0.002 


0.004 


0.007 


0.039 


Total 


Sf 


0.012 


0.017 


0.036 


0.040 


0.057 


0.087 




Cf 


0.016 


0.018 


0.020 


0.025 


0.029 


0.054 



and with the theoretical estimates of the magnitudes of decay mode in our CP sample and for the combined 
CKM matrix elements within the context of the SM [19]. J/iIjK°{K'^ + K°) mode. CP violation in rjcK^ mode 
We also report measurements of Cf and Sf for each is established at the level of 5Aa standard deviation in- 
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eluding systematic uncertainties. 
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